Anomalous Diffusion in Infinite Horizon Billiards 
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We consider tlie long time dependence for the moments of displacement {\r\'^) of infinite horizon 
billiards, given a bounded initial distribution of particles. For a variety of billiard models we find 
drl') ~ t"^' (up to factors of logt). The time exponent, 7,, is piecewise linear and equal to q/2 for 
q < 2 and q — 1 for q > 2. We discuss the lack of dependence of this result on the initial distribution 
of particles and resolve apparent discrepancies between this time dependence and a prior result. The 
lack of dependence on initial distribution follows from a remarkable scaling result that we obtain 
for the time evolution of the distribution function of the angle of a particle's velocity vector. 



I. INTRODUCTION 



Diffusion of particles in an infinite domain billiard is a well studied problem [1-5,9,14]. By a billiard we refer to the 
motion of a point particle in a two dimensional domain in which the particle moves with constant velocity in straight 
line orbits executing specular reflection (i.e., angle of incidence equals angle of reflection) from fixed boundaries. By 
an infinite domain we refer to an unbounded two dimensional region. An early consideration of diffusion in a billiard 
as a model in physics was made by Lorentz [1] to model electrons in a metal. In this model (called the Lorentz gas) 
, particles move freely and reflect specularly from fixed, randomly placed, hard-wall scatterers. A modification of the 
two dimensional Lorentz gas in which there are circular scatterers on a square lattice is an example of an infinite 
. ^ horizon billiard, called the Sinai billiard [5], and is illustrated in Fig. 0(a). Infinite horizon billiards are the subset of 
r—^ infinite domain billiards that contain channels through which a particle may pass without ever reflecting off a billiard 
I " 1 ! wall. 

In this paper we consider diffusion in infinite horizon billiards. The examples that we study numerically are shown 
^-H ■ in Figs, [^(a)-(d). The billiards in Fig. |l| include: (a) the Sinai billiard, composed of circular, hard wall scatterers 
K*" arranged on a square lattice such that the scatterers do not touch each other; (b) a modification of model (a) in which 
the circular scatterers are randomly displaced (random in direction and magnitude) by at most A < L/2 — i?, so that 
, there are channels of width L — 2{R + A) accommodating free motion; (c) randomly oriented square scatterers on a 
' square lattice; and (d) the scalloped channel, in which the domain is infinite in the y direction and bounded in the x 
I ' direction by circular arc segments, each subtending an angle less than or equal to 180°. Figure |l|(dl) shows the case 
P\| of the scalloped channel where the circular arc segments are semi-circles, while Fig. ^d2) shows the case where the 
arcs subtend an angle less than 180°. Particle motion for the situation in Fig. |l|(dl) is equivalent to particle motion 
for a stadium-type billiard [see Fig. ||(a)]; the particle motion within the scalloped channel can be folded into the 
•jH stadium billiard via reflection of the particle at a straight wall as it passes to the next cell. By cells we mean each 
(~| , portion of the scalloped channel domain between the dotted lines of Fig. |l](d) . In a similar manner particle motion in 
the bounded billiard of Fig. ^(b) can be thought of as equivalent to motion in the infinite billiard of Fig. |l|(a). 

One important means of characterizing transport in an infinite domain two dimensional billiard is through the 
phase space probability distribution function (pdf), P{x,y,9,t), where 9 is the angle of the particle velocity and we 
?-j take all particle velocities to have magnitude one. From the pdf one can calculate the displacement moments of the 
distribution. The qth moment at time t is 



{r") = {{x' + J/2)?/2) ^ r f (a;2 + y2)«/2p(^^ t)dxdy^, (1) 

Jo Jtz 

where r — \/ x^ + y^ and TZ is the (infinite) spatial domain of the billiard and where throughout this paper we take 
the initial pdf, P{x,y,6,0), to be bounded, \P{x,y,6,0)\ < K, and to be zero outside some finite region. For the 
infinite horizon billiards shown in Fig. |l] we find that the moments of the displacement have a time dependence 

(r')^r', (2) 
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FIG. 1: The four infinite horizon bilhard structures that we consider include:(a) the Sinai bilhard such that the channel 
width W — L — 2R > 0, (b) the Sinai biUiard with random displacements A away from the square matrix such that W = 
L — 2{R + A) > 0, (c) randomly oriented squares such that W — L — \/2d > 0, (d) scalloped channel with (dl) semicircular 
arcs and (d2) arcs subtending an angle less than 180°. 
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FIG. 2; (a)The stadium billiard has a domain bounded by two semi-circular arcs of radius R and two straight lines of length 
W. Particle motion in this billiard can be mapped to motion in the scalloped channel [Fig. ^(dl)]. In a similar manner particle 
motion in the bounded domain billiard shown in (b) can be mapped to motion in the infinite domain billiard of Fig. ^(a). 



which we use as shorthand for 



7? 



t^oo log t 



(3) 



For all cases in Fig. we find the exponent 7^ to be 



^ = / 9/2 g < 2 
^« \q-l q>2 



(4) 
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FIG. 3: The dependence of the gth moment, drj'), on time for the models discussed in this paper. Shown is the numerically 
determined time exponent 7,; with the symbols x, □, *, and + corresponding to the billiards depicted in Figs. 0(a), 0(b), 0(c), 
and 0(dl) respectively. The dashed line represents 7, — q/2 and the dotted line represents 7, = g — 1. 

Results of the form (||) with 7^ composed of piecewise linear functions (different from (||)) have also been obtained in 
other situations of Hamiltonian transport [6, 7]. The occurrence of an exponent 7^ 7^ q/2 is commonly referred to as 
anomalous diffusion, and 7^ > q/2 [jq < q/2) is called superdiffusion (subdiffusion). 

Figure ^ shows the results of numerical experiments testing (Q) . In these numerical experiments we start with a 
cloud of many initial conditions distributed uniformly in the accessible space occupied by one cell [the region outside 
of the scatterer and within < y < L,0 < x < L for Figs. |l|(a)-(c), and < y < i for Figs. |l|(d)] and uniformly in 
< 9 < 2tt. We then evolve the orbit of each particle in the cloud forward in time [8] and obtain (r^). In all cases 
the results we obtain are consistent with (^. The motivation for choosing an initial distribution uniform in space 
and angle is that motion of an orbit from a typical initial condition in the billiards of Fig. ^ is known to be ergodic, 
generating an invariant density uniform in space and in angle < 6 < 2tt. In another set of numerical experiments we 
used initial conditions distributed uniformly in intervals of 6 such that the lengths of the initial flights are bounded 
[e.g., for Figs. 0(a)-(c), 9o<9 <tt/2- Oq^/2 + 9^ < 9 < it ~ 9q, t: + 9q < 9 < 37r/2 - 6*0, 37r/2 + 9q<9 <2tt-9o 
where 6*0 < 7r/4]. Again, agreement with ^ and (Q) was found. This indicates that the angular particle distribution 
relaxes to the uniform distribution sufficiently fast that the results (^ and ^ are not modified. Explaining the 
reason for this insensitivity to the initial distribution is one of the main contributions of this paper (Sec. 2C). Our 
other main contribution is the result that, for long time, initial distributions with no particles in a interval about a 
direction of infinitely long flight (say 9 — Q) lead to long-time distributions with an invariant scaling form. Specifically, 
P{9,t) — J^P{x,y,9,t)dxdy approaches P{(f>), where (f) = 9t and P{(j>) is independent of t. Furthermore, the same 
P(0) universally applies for the biUiards of Figs. 0(a), 0(b), 0(c), 0(d2), but P((/)) is different for the billiard of 
Fig. 0(dl). We use this scaling result to show the insensitivity of (^) to the initial particle distribution. 

Bleher [9] showed for infinite horizon billiards that the limit in distribution as t — » 00 of the pdf of the particle 
displacements in the billiard is Gaussian with a width that increases with time as ^/tTogi. (For ordinary diffusion, 
the result is the same except that the width increases as %/I.) However, the asymptotic t dependence of the moments 
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cannot be calculated from Bleher's result. This is discussed in Sec. 3. Note that the definition of the symbol ~ 
given in Eqs. (^) and (||) is such that logarithmic corrections to the scaling of (r') with time are not included [e.g., if 
(r^) = Ctlogt, where C is a constant, as suggested by [9], then 72 from (||) is one, consistent with (^J. 



II. THEORY 



The difference between particle transport in an infinite horizon billiard and particle transport in the case of normal 
diffusion is due to the arbitrarily long trajectories found in the infinite horizon billiard. These long flights occur in the 
channels between the scattering boundaries of the billiard. When a particle is traveling nearly parallel to the axis of 
the channel [for Figs. |^(a)-(c) there are many channels parallel to both the x and y axes whereas for Fig. |l|(d) there is 
only one channel, which is parallel to the y axis], it will travel long distances between reflections off the billiard wall. 
Let 9n be the angle between the trajectory of a particle and the axis of the channel after the nth reflection of that 
particle with a billiard wall. The length r„ of the flight is of the order of 1/|6'„| for |6'„| <C 1. 

For some of the billiards we consider there are strong correlations between the 9 values from one reflection to the 
next. This is especially true of the scalloped channel with semicircular arcs (stadium billiard). Here, upon reflection, 
the angle 6'„ can change by, at most, a factor of three [12]; |0„|/3 < \On+i\ < 3|0„|. Thus if \9n\ is small, |0„+i| is 
also small, and both represent long flights. This suggests an extreme model for the particle transport in which the 
length of a flight for each particle docs not change from reflection to reflection; the length of each flight is completely 
correlated with the previous flight, but different particles have different flight lengths. 



A. The Completely Correlated Model 



Consider a one dimensional system in which an ensemble of particles executes a random walk. The particles in the 
ensemble differ from each other in the length of the step, Ar^, each particle takes. 

Ar5 = l/e, ee[0,l]. (5) 

The random variable ^ is distributed uniformly in [0, 1] is inspired by 0„ defined above, but is chosen to occupy 
the interval [0, 1] for simplicity), and ^ for a particle does not change from step to step (complete correlation). Since 
every billiard particle has the same speed we normalize the magnitude of the velocity to 1 and thus take the time 
between steps in our one dimensional random walk model to be 

A<« = 1/e (6) 

Since the random walk behaves like normal diffusion when the number of steps is large, the qth moment for each ^, 
(|r|9)j, wiU be 

(l^l'>,-(^)'/^fort»At,. (7) 



Substitution for Ar^ and At^ yields, for long time t, {t ^ At^) 



^ {t/O"^' (8) 



for each f . 

To find the qth moment for the entire ensemble one needs to take an average over all particles: 



r 



9\ 



:de (9) 



This average has two main contributions: 



{\rn- t {tliY'^Ai+ f'f^Ai (10) 

JN/t Jo 

where N is chosen to be large, so that the random walks have made many steps. The first term in Eq. (^0|) represents 
the fraction of particles that have executed at least N steps and so can be described by Eq. (||). The second term 
comes from those particles that are still in their first flight (with velocity = 1). The contribution from the interval 
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N/t > ^> 1/t [omitted in (|io|)] has neither the majority of walkers (in the hmit of long time) nor the most extreme 
displacements and so does not dominate the other terms. For long time Eq. ( [l0| ) yields 

f^/' q<2 , 
t"-! q>2 ■ ^^^^ 

This simple toy model shows that it is the very long flights allowed by the open channels that give rise to the 
non-Gaussian behavior for moments greater than 2. One might then think that for the real billiard systems of Fig. 
|l| the behavior for q > 2 is a trivial result of the initial distribution of angles. However, this is not the case and 
the behavior ( pi] ) is robust to changes in the initial distribution of particles. As discussed in Sec. 2C, non-uniform 
initial angular distributions scatter rapidly enough that ( |lll ) still holds even if the initial distribution has no particles 
traveling nearly parallel to the channels. 

B. Moment Equation for Infinite Horizon Billiards with a Uniform Initial Angular Distribution 

We now consider a uniform initial spatial (within a cell) and angular distribution and for the case of the ergodic 
billiards of Figs. |l|(a) (Sinai billiard) and ^(d) (scalloped billiard). In these cases such a distribution is stationary when 
(x, y) is taken modulo the appropriate cell period. We show that jg is given by Eq. ^ for these real billiard systems. 
The particle transport in an infinite horizon billiard must proceed at least as quickly as a random walk process (i.e., 
as fast as normal diffusion). While there exist mechanisms for faster than diffusive transport (to be discussed below), 
there is no stable mechanism to stop or trap a billiard particle. The periodic orbits that might trap the particle in 
these systems are all exponentially unstable. That a random walk is a lower bound on the transport in the billiard 
system implies that 7q > q/2. 

In addition to diffusive like behavior, particle trajectories that consist of a single long flight ("ballistic" flight) also 
participate in particle transport. If every particle were to move ballistically, we would find 7g — q. This provides 
an upper bound on 7^; not every particle will execute a single uninterrupted flight. We can, however, place a lower 
bound on the fraction of particle that do. Since the distribution of particle velocities is uniform in orientation, the 
fraction of particles executing flights of distance vt or more (where v is the velocity of the particle) involve a fraction 
of order of W/t of the particles occupying a channel at any one time, where W is the channel width defined in Fig. [|. 
(Although we have deflned v — I, we retain v in this section to clarify when we are speaking of distances and when 
we are speaking of times.) The contribution to {\r\'') from these particles is 



> 



{vty{c/t) = ct'^-\ (12) 



where C is a constant depending on the size of the channel. 

There are two points that can be fixed on the graph of 7^ vs. q. The first is the zeroth moment, which by the 
conservation of probability is identically equal to one. Thus 70 = 0. The second point is fixed by the diffusion 
coefficient, which relates the second moment to time. It has been found [14] that (r^) is of the order tlogt and so, 
consistent with the definition (||), ignoring factors of logi we have 72 = 1. 

Next we argue that 7^ is a concave up function of q. Invoking the Cauchy-Schwartz inequality (||(a;-j/)|| < ll^^ll x ||?/| |), 
for < e < q 

{\r\i) = (|r|('J+^)/2|r|(9"^)/2) < (|r|«+^)i/2(|r|«~^)i/2. (I3) 

Thus 

log(|r|''+^) + log(|r|«~^) 
log(|r|^) < ' ' ^ (14) 

From (^, (H), and ( |l^ ) we have that the graph of 7^ versus q is indeed concave up 

7. < 2i±i±^. (15) 

Finally, we show that Eq. (||) holds, i.e., that 7^ = q/2 for q < 2 and 7g = <? — 1 for q > 2. We have given lower 
bounds on 7^, i.e. 7^ > q/2 and 7g > 5 — 1- We have also pinned the value of 7^ at two values: 70 = and 72 = 1. 
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The first lower bound, the two known values of 7g, and the concavity condition (|T^ force 7g = q/2 for q <2. There 
is a trivial upper bound on 7^, 7^ < q. This upper bound, along with Eq. (|lj) implies that the slope of 7^ never 
exceeds one. Since 72 = 1, 7<; < <? — 1 for (7 > 2. Therefore, for q > 2 both the upper bound and the lower bound 
coincide resulting in = q — 1 for q > 2. 



C. Non-uniform Initial Angular Distributions 



The discussions of the previous two subsections relied on the existence of a uniform distribution for 9n in the particle 
ensemble. This applies if one starts with a distribution uniform in angle and uniform in the space within a cell, in 
which case the particle ensemble will retain the uniform angular distribution. On the other hand, even if we have an 
initial distribution that is non-uniform, it will (under very general conditions) relax to the uniform distribution as 
time increases. The question we now address is whether this relaxation is fast enough to yield the same result for the 
temporal scaling of Eq. (pT|), as the initially uniform distribution. 

To consider this question, we first study the relaxation of an initial particle distribution with no particles in a 
finite gap around the channel direction. For example, for the scalloped channel we consider the case where the initial 
particle distribution is uniform in the space within a cell C and 



P(xv9t-Q)-l^ \0\ < Oo or \9 - n\ < Oo, . 

i^{x, y, y, t - uj - I ^ otherwise ' ^^^> 

where < < 27r is the angle of the particle velocity vector with the y-axis, K and 60 are constants, and {x, y) is in 
C. Similarly, for the cases of Figs. 0(a)-|l](c) we consider an initial distribution. 



^ ' ' y K otherwise ^ ' 

In particular, we focus on the behavior of P{x,y, 9,t) for \9\ small, t large and {x,y) in a channel. (Similar results 
apply for ( |l7| ) with \9 — 7r/2| or \9 — tt\ ot\9 — 37r/2| small.) For all cases shown in Fig. |], we find a remarkable scaling 
behavior. Let P{9, t) — J J Pdxdy where the spatial integral is over a channel. Then, if we introduce the scaled 
variable ip — 9t, we find that, in all the cases we have tested, the angular distribution function P{9,t) approaches a 
stationary form. That is. 



P{9,t) ^ P((t>) a.st^ 00. (18) 

We illustrate this with a numerical calculation on the scalloped billiard with 180° arcs (Fig. [|(dl)) in Fig. ^. In 
generating this figure we follow the evolution of a large number of orbits initialized according to Eq. dlq) , and form the 
distribution P{9, t) using a histogram approximation. As shown, at successively larger t the distribution approaches 
a time independent form P{(j)) where P{(j)) is zero for (p < 4>mm, and, as (f> increases past (pmin, Pi4') increases, 
asymptoting to a constant for large (/). 

A similar numerical experiment for the Sinai billiard (Fig. |l|(a)) yields the result in Fig. In this case the large 
time distribution assumes the form. 



where C and 0o are constants. Furthermore, we numerically obtain this same form for all the other cases of Fig. |^, 
[except for the case of the scalloped billiard with 180° arcs (Fig. |l|(dl)) which gives the result in Fig. |j. We explain 
the reason for the result ( pj| ) and why it does not apply for the billiard of Fig. |l|(dl) subsequently, but before doing 
that we first show that these results for P{(f>) imply the applicability of Eq. (^ for the large time behavior of (|r|*). 

In order to show that (^ applies consider the fraction 2(io) of particles that, at some large time to, have < (/)*, 
where we take = 2(/)niin for the case of Fig. |l](dl) and 0, = (po for the other cases. Noting that </> < 0* implies 
9 < (j)^,/t, we have that for large to 



^{to) = f'' ° P{9t)A9 = {toY^ P(0)d</) = Kt-\ 
Jo Jo 



(20) 
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FIG. 4: The time evolution of an initially uniform angular distribution of particles in the interval [.3,7r/2] for the scalloped 
channel with semi-circular arcs. The distribution is shown at three times; + ai t = 27, x at 81, and * at 243. The distribution 
of particles in the channel of width W becomes self-similar and static in the coordinate (j> = 6t with P{(j>) — for (j) < (jjmin- 

where K is a constant. Between times to and 2to these particles experience flights of length ^ vto- Hence for t = 2to 
these flights give a contribution to {r'^) that is approximately [2K /t){vtY ^ t'^^. Thus the lower bound ( |l^ ) still 
applies, and, by the reasoning in Sec. 2B, we again obtain (^. 

We note that the asymptotic time dependence P{<j>) with cj) ^ 9t is marginal in the sense that, if the repopulation 
of an initially empty channel were slower (in the sense below), then (^ would not be recovered. In order to see this, 
consider the hypothetical case where an asymptotic distribution P{4>) of the form in Fig. jijor Fig. ^ is still approached, 
but with a self-similar scaling given by (pa ~ Ot"^ . We have already considered the case a = 1, while a < 1 (a > 1) 
corresponds to slower (faster) filling in of the channel. If, for < a < 1 (i.e. slow repopulation) we pursue the same 
reasoning as above for the a — 1 case, then we obtain the bound jq > a{q — 1), and we can no longer conclude that 
(^ holds. For a > 1, we consider at time to a range of angles 6^ 1/to > > 6^ ^ 'f>*lto- Since a > 1, we have 
that 6^ 9^, and replacing 0_ by zero does not alter the estimate for the contribution to (?■'') from 9+ > 9 > 9^. 
Thus the lower bound estimate of Sec. 2B still applies for < a < 1. 

We now discuss the asymptotic forms P{<p) illustrated in Fig. ^ for the scalloped billiard with 180° arcs, and in 
Fig. H and Eq. ( pj| ) for the other cases. 

First we discuss the scalloped billiard with 180° arcs. A full treatment of the theory yielding P{4>) in this case 
will be given elsewhere [13]; in the present paper, we will limit our discussion to the basic reason for the difference 
between this case and the other cases. In particular, we discuss why, for this case, P(0) = in a finite interval about 
zero, (p < 0inin- From Fig. |l|(dl) we see that after a long flight in the channel a particle will collide with the channel 
wall close to one of the cusp points where two arcs touch. For the case of 180° arcs, the tangent to such a section of 
the channel wall is nearly horizontal. Thus, upon reflection 9 = min(|0|, \9 — 7r|) will still be small. In fact as shown 
in [12] by consideration of the geometry, ()n+i on the (n -|- l)st reflection cannot change by more than a factor of 3 
from 9r, , 



On/3 < On+i < 34, for 0„ < 1. 



(21) 



Thus, if 9n is small, 9n+i is still relatively small. We can obtain the lower bound on (p by considering the most extreme 
case where 9 always decreases by 3 on every bounce. 



(22) 
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FIG. 5: The time evolution of an initially uniform angular distribution of particles in the 6 interval [.3, 7r/2] for the Sinai billiard 
with L/R = 1. The distribution is shown at four times; + at f = 9, x at 27, * at 81, and □ at 243. The distribution of particles 
in a channel becomes self-similar and static in the coordinate A — 9t. 



W/{v9n), (23) 



where W is the channel width defined in Fig. |^(dl), v is the particle velocity, and i„ is the time of the nth reflection. 
Multiplying (^) by (p3|) we have 



$n+l = ^4>n + W/V, (24) 

which, for large n, asymptotes to the solution (j> — 3W/(2v). Thus (pnun = 3W/{2v), which agrees with our numerical 
solution Fig. ^ (see also [13]). 

In contrast to the case of the scalloped channel with 180° arcs, in the other cases shown in Fig. the scattering 
of a long flight upon reflection from a channel wall leads to a much more drastic change in the angle of a particle's 
velocity vector with respect to the channel axis. For example, for the case of the Sinai billiard, the angular deflection 
is typically of order 9^^^ which , for small 9, is much larger than 9. For the case of the scalloped channel with arcs of 
less than 180°, a particle moving nearly parallel to a channel axis is scattered by an angle of order one. Furthermore, 
after a large deflection, the orientation of the particle's velocity vector is rapidly randomized by a succession of 
many reflections which, since the particle is no longer in a long flight, can occur in a relatively short time. These 
considerations lead us to a model for the cases in Figs. 0(a)-|l|(c) and |l|(d2) in which we adopt the model hypothesis 
that, when a particle in a long flight suffers a collision with a billiard wall, the orientation of its velocity vector is 
randomly scattered with uniform probability density in [0, 27r]. We wish to determine the evolution from the initial 
condition (a) in Fig. ^for the case f?max ^ 1- This initial distribution is equal to the initial distribution (b) in Fig. ^ 
minus the initial distribution (c) in Fig. ^. The initial condition (b), which is uniformly distributed in angle, remains 
unchanged when it is evolved forward in time [Pi,{9,t) — Pi,{9,0)], since it is an invariant distribution. Thus, to find 
the evolution from initial condition (a), we can determine the evolution from (c), and then subtract it from (b). The 
long time evolution from (c) can be found by considering the time at which particles are scattered. Consider, for 
example, the scalloped channel. Fig. |l|(d2), and a particle with a small initial 9o- Suppose the particle is located in 
the channel at a distance Ax from the boundary of the channel with which it will collide [left or right vertical dashed 
line in Fig. |l|(d2)]. If Ax < vt sin 9o = vt9o, the particle scatters; if Ax > vt9o, it does not scatter. Since the particles 
we are considering are in the channel. Ax < W, every particle with 9o > W/ (vt) must have scattered at least once. 
We assume that t > W / {v9max) = ^o- Since 0max is small, the scattered particles contribute a small positive value 
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FIG. 6: An initial distribution of particles with trajectories uniform in the angle 9 for 9 > ^max and zero for angles 9 < 6'max, 
Fig. (a), is equivalent to a distribution uniform for all 9, Fig. (b), minus a distribution that is uniform for 9 < 6'max and zero 
for 9 > Smax, Fig. (c). We use the time evolution of the distributions in Figs, (b) and (c) to find the time evolution of Fig. (a), 
see Fig. ^. 




of order 6'max to PdO^t) in < 6* < 27r. Thus Pci0,t) is small (i.e., of order 6'max) for 9 > W/vt. For 9o < Ax/(ut), 
t > to, the particle has not yet scattered. Assuming that the initial spatial distribution of particles in the channel is 
uniform, the fraction of particles with initial angle 6o that have scattered is doVt/W. Thus 

where we have neglected the small, order 9max, contribution to Pc{0,t) from scattered particles. Subtracting Pc from 
Pb as illustrated in Fig. ^ we obtain the time asymptotic form in Fig. ^ and Eq. (p^. 
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III. DISCUSSION 



As mentioned in the introduction, Bleher [9] proved, up to some natural conjectures that "For any periodic config- 
uration of scatterers with an infinite horizon the limit in distribution 

t^oo (t log 1/2 ' ^ ' 

exists and 77 is a Gaussian random variable." (We have substituted our r{t) for Bleher's x(t) for the sake of notational 
consistency) . 

We take no issue with Bleher's proof referred to above, but we emphasize its use of limit in distribution convergence. 
This kind of convergence is the least restrictive kind of convergence considered within the context of probability and 
statistics. A proof of convergence in distribution implies only the ability to calculate the expectation value of functions 
that remain bounded. The convergence required by a limit in distribution is only strong enough to allow calculation 
of the expectations of functions that remain bounded [10]. Therefore the convergence of Bleher's pdf is not strong 
enough to allow the moments of the distribution to be calculated. As a result of the strong weight \r\'^ puts on the tails 
of the distribution (|r| — )■ 00), two different distributions with the same limit in distribution can have very different 
moments. 

The fact that Bleher is able to accurately (according to our simulations) calculate the second and lower moments 
of the displacement distribution suggests that his result can be strengthened to "convergence in q*^ mean" which is 
satisfied for a sequence Xn if the expectation value of |X„ — ATj'^^Oasn— >oo. Convergence in q*^ mean also implies 
that the expectation value of |Ar„|P limits to the expectation value of \XY' for \ < p < q [15]. Thus our results are 
consistent with convergence in gth mean to Bleher's distribution for 9 = 2, but rule out convergence for any higher 
value of q. 

The inapplicability of Bleher's result explains the discrepancy between Eqs. (|^)-(|^) and the result {\r\'^) ^ (i logt)'/^ 
one would find by mistakenly calculating moments using Bleher's pdf. Equations (|2|)-(||) also differ from normal 
diffusion, (|r|9) - i'/^ as weU as fr om the result suggested in [11]. 

In conclusion, our two main results are as follows: 

(a) The moments (r') scale as t'^'' with 7^ given by Eq. (^) for any initial bounded distribution, |P(a;, y, 6*, 0)| < K, 
that is zero outside some finite region (in particular, Eq. (^ still applies if the initial distribution has no particles 
with infinite flights). 

(b) If the initial distribution has no particles in a 0— interval about a direction of infinitely long flight (say = 0), then 
P{0,t) — J^P{x,y,9,t) approaches a time-invariant scaling form P{4>), where (f> = 9t and P{4>) is universally 
the same (Fig. |) for the biUiards of Figs. |(a),|l|(b), 0(c), |l|(d2), but is different (Fig. |) for the bilhard of 
Fig. |(dl). 
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